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^ ! Abstract 

By employing the perturbative QCD (pQCD) factorization approach, we calculated some 
important next-to-leading-order(NLO) contributions to the two-body charmless hadronic decays 
B^ — > p'^Tj^'' and B^ — > p^ {lo , (j))ri^'\ induced by the vertex QCD corrections, the quark-loops 
as well as the chromo-magnetic penguins. From the numerical results and phenomenological 
r~| . analysis we find that (a) for B — > p rj^'' decays, the partial NLO contributions to branching 

♦^1 ratios are small in magnitude; (b) for B^ -^ p^ {uj , <j))ri^'' decays, the NLO contributions can 

Qh| provide significant enhancements to the leading order predictions of their branching ratios; and 

(— I ■ (c) the pQCD predictions for the CP-violating asymmetries A'(j'^p{B^ — > p^rj^'^) are consistent 

with the data, while the predicted Acp{B^ -^ p^{uj)r]^'>) are generally large in magnitude and 
^SJ . could be tested by the forthcoming LHCb experiments. 
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I. INTRODUCTION 

During the past decade, the B factory experiments have achieved great successes. More 
than one bilhon events of BB production and decays have been accumulated and analyzed 
by BaBar and Belle Collaborations. The forthcoming LHC experiments will provide 2-3 
orders more B meson events than the B factory, and high precision measurements for the 
branching ratios and CP- violating asymmetries of many B meson rare decays will become 
true within the following three to five years. Now it becomes a very important and urgent 
task to reduce the uncertainty of the theoretical predictions, in order to test the standard 
mode (SM) and to find signals or evidence of the new physics beyond the SM through 
the B meson experiments [jj]. 

For the charmless decays B — >■ M1M2 ( here Mi are light mesons composed of the 
light u,d,s quarks), the dominant theoretical error comes from the large uncertainty in 
evaluating the hadronic matrix elements {MiM2\0i\B) . In order to increase the accuracy 
of the SM predictions, various factorization approaches have been proposed in recent years. 
The perturbative QCD (pQCD) factorization approach |3] , together with the so-called 
QCD Factorization (QCDF) [3| and the SCET [^ , are the most popular factorization 
approaches [3, ^ being used currently to calculate the hadronic matrix elements [3|, |5|, la, 

mMMMMMMM- 

When compared with the QCDF or SCET factorization approaches, the pQCD ap- 
proach has the following three special features: (a) since the kr factorization is employed 
here, the resultant Sudakov factor as well as the threshold resummation can enable us to 
regulate the end-point singularities effectively; (b) the form factors for B —>■ M transition 
can be calculated perturbatively, although some controversies still exist about this point; 
and (c) the annihilation diagrams are calculable and play an important role in producing 
CP violation. 

Up to now, almost all two-body charmless B/Bg — ^ Mi Mo decays have been calculated 
by using the pQCD approach at the leading order (LO) p, |7|, |8|, |9|, [lO|, lul, 112, |l3|, \M ■ 



Very recently, some next-to-leading (NLO) contributions to B ^ Ktt and several B 



PV decay modes [15|, [IGj have been calculated, where the Wilson coefficients at NLO 
accuracy are used, and the contributions from the vertex corrections, the quark loops and 
the chromo-magnetic penguin operator O^g have been taken into account. As generally 
expected, the inclusion of NLO contributions should improve the reliability of the pQCD 
predictions. 



In previous papers [10|, lUl], the authors calculated the branching ratios and CP violat- 
ing asymmetries of the B -^ p(uj, 0)?7*^') decays by employing the pQCD approach at the 
leading order. Following the procedure of Ref. [l5|, we here would like to calculate the 
NLO contributions to the B — > p(uj, (p)!]^'' decays by employing the low energy effective 
Hamiltonian and the pQCD approach. 

The remainder of the paper is organized as follows. In Sec. II, we give a brief discussion 
about pQCD factorization approach. In Sec. Ill, we calculate analytically the relevant 
Feynman diagrams and present the various decay amplitudes for the studied decay modes 
in leading-order. In Sec. IV, the NLO contributions from the vertex corrections, the 
quark loops and the chromo-magnetic penguin amplitudes are evaluated. We show the 
numerical results for the branching ratios and CP asymmetries of i? — > p{uj, (p)!]^'^ decays 
in Sec. V. The summary and some discussions are included in the final section. 



II. THEORETICAL FRAMEWORK 

Based on the pQCD factorization approach j^], the decay amphtude A{B — > M1M2) 
can be written conceptually as the convolution, 

A{B ^ MiM2) r^ fd%d%d% TT[CmB{ki)^MAk2)^M,{h)H{kuk2,h,t)], (1) 

where fcj's are momenta of light quarks included in each meson, and Tr denotes the 
trace over Dirac and color indices. C{t) is the Wilson coefficient evaluated at scale t. 
The function H{ki,k2,k3,t) describes the four quark operator and the spectator quark 
connected by a hard gluon and could be calculated perturbatively. The function $s and 
$Mi are the wave functions of the initial heavy B meson and the final light meson Mj, 
which describe the hadronization of the quark and anti-quark into the mesons. While the 
hard kernel H depends on the processes considered, the wave functions $5 and ^^i are 
independent of the specific processes. 

In the B meson rest-frame, it is convenient to use light-cone coordinate {p'^,p~, Pt) 
to describe the meson's momenta: p^ = -yn{p^ ± p^) and px = (p\p^) • Using these 
coordinates the B meson and the two final state meson momenta can be written as 

Mb, ^ Mb , n , Mb, o , 

Pb = -|(1,1,0t), Py = -|(1,4,0T), Pp = -|(0,l- 4, Ot), (2) 

respectively, here ry = uiv/Mb with V = p,uj or (p. The light meson (P = r]^'^) mass has 
been neglected. For the B — >■ VP decays considered here, only the vector meson's longitu- 
dinal part contributes to the decays, and its polarization vector is e^ = ^^ (1, —ry, Ot). 
Putting the anti-quark momenta in B, V and P mesons as ki, /c2, and k^, respectively, 
we can choose 

fci = (xiP+,0,kiT), A;2 = (a;2P2+,0,k2T), ^3 = (0,a;3P3", kax). (3) 

Then, the integration over k^ , /c^, and k^ in eq.([T]) will lead to 

A(B ^ PV) ~ ld.,i...d.,t,dt,t,dt,Mb. 

■Tl' [C(t)'S>B(xi,bi)'S>v{x2,b2)9p(x3,b3)H(x„b„t)S,(xt)e-^<'>] , (4) 

where bi is the conjugate space coordinate of kix- The large logarithms (Inmw/t) com- 
ing from QCD radiative corrections to four quark operators are included in the Wilson 
coefficients C{t). The large double logarithms (In^Xj) on the longitudinal direction are 
summed by the threshold resummation 17|, and they lead to St{xi) which smears the 



J3 



end-point singularities on Xj. The last term, e"'^*-*-*, is the Sudakov form factor which 
suppresses the soft dynamics effectively |2(]. 

For the studied B -^ Vr]^'^ decays, the weak effective Hamiltonian -ffe// can be written 
as fl? 



q=u,c 
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(Ci(/i)0f (^) + C2{^i)0l{^,)) + > ^ q(^) o,(/x) 



i=3 



(5) 



where Gp = 1.16639 x lO^^GeV^^ is the Fermi constant, and Vij is the CKM matrix 
element, Ci{fi) are the Wilson coefficients evaluated at the renormalization scale fi and 
Oi{n) are the four-fermion operators. For the case of 6 — >■ s transition, simply make a 
replacement of (i by s in Eq. (j5]) and in the expressions of Oi{^) operators, which can be 



found easily for example in Refs. |12l . Il8 



In PQCD approach, the energy scale "t" is chosen as the largest energy scale in the 
hard kernel H{xi, bi, t) of a given Feynman diagram, in order to suppress the higher order 
corrections and improve the reliability of the perturbative calculation. Here, the scale "t" 
may be larger or smaller than the ruf, scale. In the range of t < rrih or t > rrih, the number 
of active quarks is Nj = 4 or Nf = 5, respectively. For the Wilson coefficients Ci{fi) and 
their renormalization group (RG) running, they are known at NLO level currently |18| . 
The explicit expressions of the LO and NLO Ci{mw) can be found easily, for example, in 
Refs. 0,[l3- 



When the pQCD approach at leading-order are employed, the leading order Wilson 
coefficients Ci{mw), the leading order RG evolution matrix U{t, m)^'^^ from the high scale 
m down to t < m ( for details see Eq. (3.94) in Ref. [M]), and the leading order asit) are 
used: 

4:71 

where po = (33 - 2Nf)/3, AgJ,^ = 0.225Ge1/ and AgJ,^ = 0.287 GeV. 

When the NLO contributions are taken into account, however, the NLO Wilson coef- 
ficients Ci{mw), the NLO RG evolution matrix U{t,m,a) (for details see Eq. (7.22) in 



Ref. [18|) and the as{t) at two- loop level are used: 



where po = (33 - 2iV/)/3, /3i = (306 - 38A^/)/3, Ag^^ = 0.225 GeV and Ag^^ = 0.326 
GeV. 



III. DECAY AMPLITUDES AT LEADING ORDER 

In the pQCD approach, the Feynman diagrams as shown in Fig. [T]may contribute to 
B — >• ofg;, (t))r]^'^ decays at leading order. These decays have been studied previously in 
Refs. 10, lul by using the pQCD approach. In this paper, we focus on the calculations of 



some NLO contributions to these decays in the pQCD approach. We firstly recalculated 
and confirmed the previous calculation. For the sake of completeness, we present the 
relevant LO decay amplitudes in this section. 

At the leading order, the total decay amplitudes for B -^ prj, B^ — > ur], and B^ -^ (prj 
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FIG. 1: Feynman diagrams which may contribute to the B 
order. 
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(10) 



M{(t)r]) = -FJ^^t f as + as - i^aj - ^0.9 j ^i(0) 



- (Ma + Ma4,) 6 (C^ - Ic^o) ^2(0) 

- (Mf + Mf/) 6 (a - iCg) F2(0), (11) 

where ^„ = V*fjVud, ^t = Vtl^td, and -Fi(0), -^2(0) are the mixing factors oi r] — rj' system 
as define in Eq. l\A9\i : 

Fi(0) = — COS0, F2(0) = -sin0. (12) 

V2 

The Wilson coefficients a, appeared in the expressions of the total decay amplitude are 
the combinations of the ordinary Wilson coefficients Ci{fi), 

ai(/i) = C2{fJ,) + -Ciifi), 02 (/i) = Ci(/i) + -C2(/i), 

aiifi) = Q{fi) + ^^^^^Jor z = 3,5,7,9; or 4,6,8,10. (13) 



The individual decay amplitudes Fep,--- , as given in Eqs. flHl fTTj) . are obtained by 
evaluating individual Feynman diagrams for a given decay mode and can be written as 

Fev = Av2iiG pC pfn^B / dxidx2 / 6i(i6i62'^&2 0_B(a^i, &i) 

Jo Jo 

X { [(1 + X2)0y(x2) - (1 - 2a;2)ry(0y(x2) - 0U^2))] Eeita)heiXi,X2, h, 62) 

-2ry0t'fe)^e(C)/ie(a:2, a;i, 62, h)} , (14) 

F^ = 8V2GF'KCpm'^ / dxidx2 / hidhih2dh2(l)B{xi) 

Jo Jo 

X {-^77 [0y(^2) - ry((2 + X2)(l)v{x2) + a;20*y(x2))] Ee{ta)K{Xi, X2, bi, 62) 

+2rvrr^4>v{x2)Ee(t'^)he{x2, Xi, 62, &l)} , (15) 

Mev = M[y = —=GFTTCFm% / dxidx2dx3 / bidbibsdbs(f)B{xi,bi)(j)^{x3) 
v3 JO Jo 

X {[-X20v(^2) +2X2ry0*^(x2)] -Eg(t5)/l„(xi,X2, X3, 61,63) } , (16) 

May = —=GF7rCFm'^ / dxidx2dx'i / bidbib^db^cpsixijbi) 
v3 Jo Jo 

X { [(1 - X2)0y(x2)0;^(x3) + ryr^(l - X2) (0yfe) + 0y(^2)) (^^(^s) " ^^'(^s)) 
+rvrnX3 (0y(^2) - 0*y(^2)) (^^'(^s) + ^^(^s))] E'^{tc)hna{Xl,X2, X3, 61, 63) 

- [x30y(x2)</);J(%) + 4ryr^0f.(x2)0j'(^3) " '"y'"r,(l " a^s) (0y(^2) + 0*y(^2)) 

• (0j'(^3) - (pJiixs)) - ryr^X2 (0y(^2) - 0*y(^2)) (^^'(^3) + (pj^ix^))] 
xKit'c)Kai^l,^2,X,,b,,bs)} , (17) 



1 f\ r^ /*oo 

M^Y = —=GF'^CFm'^ / dxidx2dx3, / 61^6163^63 0B(a;i, 61 
v3 Jo Jo 



V ^ JU J u 

X { [(1 - X2)ry0;J(x3) (0y(x2) + ^'yfe)) " ^ry^^s^yfe) (^^'(^s) " ^^(^s))] 
xE^(tc)/i„a(a:i,X2,X3,6i,63) - [-(^s + I)ry0;^(x3) (0y(x2) + 0*y(x2)) 

r^(a;3- 2)0^(^2) (^^'(^s) - ^J'l^s))] ^l(t'J/i'„„(a;i,a;2,a;3, 61, 63)} , ( 



1^ 



1 fi z*! /•oo 

M^ = —^Gp'^CFm'^ / dxidx2dx-i I 61^6163^63 ^^(a^i, 61) {[(x2 — 1) 
v3 Jo io 

X0y(x2)0;^(x3) - 4r^ry0^(x2)0j'(x3) + r^ryO^s (^yfe) + ^yfe)) 

■ (0^(a;3) - 0^(X3)) + ryr,(l - 3:3) (0^.(2:2) - M^^)) ^(0:3) + ^^^(a^s))] 

■E'^(te)hna{Xi, X2, X3,bi, bs) + [a;30y (:r2)0;^(x3) +X3ryr^ (0y(^2) +0y(S2)) 

• W(^3) - 0j^(^3)) + ryr^(l - 2:2) (0y(^2) - 0U^2)) (^^(^s) + ^^(^s))] 
xK(t'J/i;,(xi,a;2,a;3,6i,63)}, (19) 

/•l /"OO 

Fe = 4v2G'i?7rCi7m'| / dxidx2 / 61^6162^62 0_b(xi, 61) 

Jo Jo 

X { [(1 + a;2)0;^(:r2) + (1 - 2a;2)r,,(0j'(^2) - 02^(^2))] Ee{ta)he{xi,X2, 61, 62) 
+2r^0j'(a;2)Fe(t'J/ie(a:2,xi,62,6i)} , (20) 

Me = ^pG'iT'TrCi???!^ / dXidX2dX3 / 61^6163^63 0b(Xi, 6i)0y(x3) 

v3 Jo io 

X {-a;20;^(^2) -2a;2r^0^(S2)} ■ -Eg(tfe)/i„(xi, a;2,X3, 61, 63), (21) 

1 28 r^ r°^ 

M^^ = -^TrCprn^^ry / dxidx2dx-i I 61^6163^63 ^^(a^i, &i) 

V 6 Jo Jo 

X {(1 - X3)0;J(X2) ■ (0y(^3) - 0y(^3)) 

+r^(l - x^) {(t)^{x2) + 0^(2:2)) (0y(a:3) - ^U^a)) 

+r^a;2 (0j'(^2) -05^(^2)) (0y(^3) +0y(^3))} ■Fe(t6)/i„(a:i,X2,X3,6i,63), (22) 

1 f\ r^ /'OO 

Ma = ^pv2Gi?7rCi?m^ / dxidx2dxz / 61^6163^63 ^^(xi, 61) 
v3 Jo io 

., r r/1 ™ \j.A/- \ i /- \ », », /I _ \ (iPi^ \ , iT/- 



— —\/Z,(jpli(_ypllLj^ I UXIUX2UX3 I flUUlf3UU3 (/^iJ^^Xl, flj 

v3 JO Jo 

X { [(1 - X2)0;^(x2)0y(x3) - r^ry(l - X2) {(t)^{x2) + 0j'(^2)) (0y(^3) " 0y(^3)) 
-rr,rvX3 (^^(^2) - 0j'(^2)) (^yfe) + 0U^3))] E'^{tc)Ka{,Xi,X2, X3, 61, 63) 

- [a;30;^(x2)0y(x3) - 4r^ry0j'(S2)0y(^3) + r,,ry(l - X3) (^^'(^2) + (fir,{x2)) 

■ {<Pv{xz) - <Pv{x^)) + r^rvX2 {<P^{,X2) - ^^(Xs)) (0y(:C3) + 0y(x3))] 

x^:(OC(a^i, 2:2, 0:3,61, 63)} , (23) 



1 ^ z*! /"Oo 

M^^ = —^GFT^CFrrig / dxidx2dx^ I 61^6163^63 0b (xi, 61) 
v3 Jo io 

X { [-(1 - X2)r^0y(x3) (0J'(X2) + 0J(X2)) - ryX30;J(x2) (^ylxg) - (pv^X^))] 
xE'^{QhnaiXuX2, X3, 61, 63) - [(^2 + l)r^0y(:r3) (0j'(x2) + 05'(:^2)) 

-ry(x3 - 2)0;^(a;2) {(pU^s) - M^s))] K{Qh'^^{xi,X2,x^MM)] , (24) 

M^^ = —j=G pT^C Frrv'B / dxidx2dx3 / 61^6163^63 0s(xi, 61) {[(x2 — 1) 
v3 Jo io 

X0;J(x2)0y(x3) + 4ryr^0j'(x2)0y(x3) - ryr^X2 (0j'(x2) + 0j'(S2)) 

■ (0t'(x3) - (fvixs)) - r^ry(l - X3) (^^(sa) - ^J'fe)) (^yl^s) + ^M^s))] 

■E'^{te)hna{Xi, X2, X3, 61, 63) + [x30;^(x2)0y (X3) - X3r^ry {(p^ {X2) + 0^(X2)) 

■ {(Pvixs) - 0y(^3)) - r^ry(l - X2) {4>^{X2) - 0j(^2)) (^yl^s) + ^U^s))] 

xE:(t;)/i;,(a;i,X2,X3,6i,63)}, (25) 

Here ry = mv/mB is the mass ratio with {mv = rup, m^, m^)] Cp = 4/3 is a color factor. 
The evolution functions E{ti) and hard function hj{xi,bi) are displayed in Appendix [Bl 

The decay amplitudes for B -^ prj', B^ -^ urj' , and B^ -^ (pi]' decays can be obtained 
easily from Eqs.® to (fTTl) by the following replacements 

Fi(0) — . F{(0) = i=sin0, (27) 

F2(0) -^ F^{<P) = cos (P. (28) 

IV. NEXT-TO-LEADING CONTRIBUTIONS 



The power counting in the pQCD approach [15| is different from that in the QCDF 
approach|3|. Here the term NLO means that the decay amplitude is proportional to 
a'^ifi). We here indeed consider the partial NLO contributions only: those from the vertex 
corrections, the quark-loops and chromo-magnetic penguins. The NLO contributions from 
hard-spectator and annihilation diagrams are not known at present. When compared 
with the previous LO calculations in pQCD [l^, the following NLO contributions will be 
included: 

1. The LO Wilson coefficients Ci{mw) will be replaced by those at NLO level in NDR 
scheme [l8|. As mentioned in last section, the strong coupling constant as{t) at 
two-loop level as given in Eq. ([7]), and the NLO RG evolution matrix U{t, m, a), as 



defined in Ref. 18 , will be used here: 



U(mi, 1712, «) = U(mi,m2) + -—R(mi,m2) (29) 

Air 

where the function U{mi, 1712) and R{mi, 7712) represent the QCD and QED evolution 
and have been defined in Eq. (6.24) and (7.22) in Ref. [l8[. We also introduce a 
cut-off Ho = 1 GeV for the hard scale "t" in the final integration. 
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FIG. 2: NLO vertex corrections to the factorizable amplitudes. 



A. Vertex corrections 



The vertex corrections to the factorizable emission diagrams, as illustrated by Fig. [21 
have been calculated years ago in the QCD factorization appeoachjsl, [l9|]. According to 
Ref . 15| , the difference of the calculations induced by considering or not considering the 
parton transverse momentum is rather small, say less than 10%, and therefore can be 
neglected. Consequently, one can use the vertex corrections as given in Ref. jl9| directly. 
The vertex corrections can be absorbed into the re-definition of the Wilson coefficients 
ai{n) by adding a vertex-function Vi{M) to them (sl. llQl 



a,{fi) -> a,(/i) + ^C^^\/,(M), for ^ = 1,2; 



ttjifi) -^ aj{fi) + 



Air 3 



An 



-Cf 



N. 



V,{M), for J 



10, 



(30) 



where M is the meson emitted from the weak vertex. When M is a pseudo-scalar meson, 
the vertex functions Vi{M) are given ( in the NDR scheme) in Refs. [l5|, Il9l |: 



V,{M) 




2V2Nr 



IodxK{x)g{x) 



18+ , 

JM 



for 



2y/2 Nc 

h 



Jm ^ 



Jq dx(f)^,j{x)h{x) 



for 
for 



1-4,9,10, 

2 = 5,7, (31) 

i = 6,8, 



where /m is the decay constant of the meson M; 4>m{x) and 0^/(a;) are the twist-2 and 
twist-3 distribution amplitude of the meson M, respectively. For a vector meson V, 



^m{4>m) is replaced by (pvi.'Pv) and /m by fy in the third line of the above formulas. 



The hard-scattering functions g{x) and h{x) in Eq. flHTl) are: 



9[x} 



1 -2x 



1 



In X — in 



X 



+ 



2 In T* 

2Li2(x) — In^ x H (3 + 2iTT) \nx — {x ^^ 1 — x) 

1 ~ X 



h{x) = 2Li2{x) — In^ X — (1 + 2iTc) lux — {x ^^ 1 — x) 



(32) 
(33) 



where Li2{x) is the dilogarithm function. As shown in Ref. [I5j, the /^-dependence of 
the Wilson coefficients aj(yu) will be improved generally by the inclusion of the vertex 
corrections. 
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FIG. 3: Quark-loop diagrams contributing to i? ^ p{Lo)ri^'' decays. 



B. Quark loops 



The contribution from the so-called "quark-loops" is a kind of penguin correction with 
the four quark operators insertion, as illustrated by Fig. [31 In fact this is generally called 
BSS mechanism [20], which plays a very important role in producing the CP violation in 
the QCDF/SCET approaches. We here include quark-loop amplitude from the operators 
Oi2 and 03_6 only. The quark loops from Oj^iq will be neglected due to their smallness. 

For the b ^ d transition, the contributions from the various quark loops are described 
by the effective Hamiltonian H^ll 15| . 

Crp 



H 



ill) 

eff 



E E ^^.^^^ 



q=u,c,t ql 



:,^ ^'(^^(/i, I') (rf7p (1 - 75) T%) iWTW) , (34) 



where P being the invariant mass of the gluon, which connects the quark loops with the 
q'q pair as shown in Fig. [3l The functions C^'^\jj, P) can be written as 



C<^)(^,/' 



G^'\f^j')-I 



C2(/i), 



(35) 



for q = u,c and 



G^'\fi,f 



Csifi)+ Yl G^'^"\f^,f)[C,{fi) + CM]. (36) 

q//=u,d,s,c 

The integration function G^''\fi, P) for the loop of the quarks q = (m, d, s, c) is defined as 



G'^''\^,P 



A I dx x{l — x) In ■ 



m„ 



x{l — x)P 



(37) 



'0 ^^ 

where ruq is the quark mass. The explicit expressions of the function G^'^^fx, P) after the 
integration can be found, for example, in Ref. [ISj]. 

It is straightforward to calculate the decay amplitude for Fig|3t and [3)d. For the case 
oi B ^ V OT B ^ rj transition, we find two kinds of topological decay amplitudes: 



M, 



(9) 



Vri 



-16m^- 



C2 



dxidx2dx3 / bidbib2db2(f>B{xi,bi) {[{1 + X2)(l>v{x2)(t)^{x3] 



-rv (1 - 2x2) {(Pv{x2) - 4>v{^2)) (p^ixs) - 2r^M^2)<P^{x3) + 2rvr^ ((2 + X2) 

■M^2) + X2(p'vix2)) <(X3)] E^'^\tq,f)K{x2,Xi,b2,bi) + [-2ry0^^(x2X(x3) 

+4ryr,0^(x2)0j'(x3)] E'^'^\t'^,l'^)K{xi,X2,b,,b2)} , (3^ 
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for B ^ V transition, and 

M^y = ■-=GFCpm'^Q / dxidx2dx3 / bidbib2db2 4>B{xi,bi) [[{1 + X2) (p^{x2)4>v{x3) 

v3 Jo Jo 

+r^ (1 - 2x2) {(p^{x2) - 0j'(^2)) 0v(^3) - 2rv0;^(x2)0y(x3) - 2r^rv ((2 + X2) 

■(^^{X2) + X2(t)^{x2)) 4>v{x3)] E^''\tgJ^)he{x2,Xi,b2,bi) + [2r^0j'(x2)0y (xg) 

-4r,rv^0j'(x2)0t^(x3)] i^^-^H^;, ^")/^e(xi,a;2, 61, 62)} , (39) 

for B ^ 7] transition. Here V represents p, ui, or meson, and r^ = nij^/mB-, ry = my/mB- 
The evolution factors take the form of 

E('?)(t,/2) = C^''\t,P)al{t)-exp[-Sat], (40) 

with the Sudakov factor Sab and the hard function he{xi,X2,bi,b2) as given in Eq. (1B9P 
and Eq. (JB2I) respectively, and finally the hard scales and the gluon invariant masses are 

tg = max{y/x^mB, ^/x^mB, ^/iT^^X2Jx^mB,l/bi,l/b2); , 

t'g = max{y/x^mB,^/a^mB,^/\x^^^x^\mB,l/bl,l/b2), (41) 

/2 = (1 _ X2)x3m% - |k2T - ksxl^ ^ (1 - X2)x3m%, 

l'"^ = {x3 - xi)m| - |kiT - ksxl^ ^ (xs - xi)m|. (42) 

For i? ^ V?7' decays, we find the similar results by making appropriate replacements, 
such as r^ — > r — t]', etc. 

Finally, the total "quark-loop" contribution to the considered B -^ Vr]^''^ decays with 
V = p,uj can be written as 

Mi^^i, = < Vv^'^mi^lB >= Y: K [m^^Io+M;^1] , (43) 

q=u,c,t 

where Xg = VgbV*^. The quark-loops do not contribute to -B — » (f)-!]^''^ decays. 

It is note that the quark- loop corrections are mode dependent. The assumption of a 
constant gloun invariant mass in FA introduces a large theoretical uncertainty as making 
predictions. In the PQCD approach, the gluon invariant mass is related to the parton 
momenta unambiguously. 

C. Magnetic penguins 

As illustrated by Fig. HI the chromo-magnetic penguin operator Ogg also contribute to 
B -^ Vr]^'' decays at NLO level. The corresponding weak effective Hamiltonian contains 
the b ^ dg transition, 

K7/ = -^V,bV:, Cf/ Os,, (44) 

with the chromo-magnetic penguin operator, 

Os, = ^rrib ~d,a^^{l + l,)Ttp%b,, (45) 
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B 



y(?7W) B 



(a) (b) 



FIG. 4: Chromo-Magnetic penguin (Osg) diagrams contributing io B ^ p{Lo)ri^'' decays, 
where i,j being the color indices of quarks. The corresponding effective Wilson coefficient 



Cgg — Csg +C5 [15 



In Ref. [2IJ, the authors calculated the chronio-magnetic penguin contributions to 
B ^ (f)K decays using the pQCD approach. They considered nine chromo-niagnetic 
penguin diagrams corresponding to the non-local operator Og^, as given in Eq. (2.3) of 
Ref. [2l|, generated by operator Ogg as defined in Eq. (HSj) . The first two Feynman 
diagrarns (a) and (b) in Ref. |2l(] are the same as Figs. UK and Hb here. According to 
Ref. \21l\, the diagrams (a) and (b) dominate, while other seven diagrams are small or 
negligible. It is therefore reasonable for us to consider the NLO contributions induced by 
the diagrams (a) and (b) only, for the sake of simplicity. 

The decay amplitude for Figs. SK and Sb can be written as 

M\^'^ = 16m|— ^= / dxidx2 dx^ / bidbib2db2 0b(xi, 61) {[- (1 - X2) {20^(^2) - ry 
2v A'c Jo Jo 

■(30y(x2) - (t)v{x2)) - rvX2 {(t>v{x2) + (t)v{x2)) } (t>^{x3) + r^ (1 + X2) X3(t)v{x2) 

■ (3<(X3) + <(X3)) - ryr, (1 - x^) {<PU^2) + 0M^2)) (3<(X3) - <(a;3)) 

-ryr^Xs (1 - 2X2) (0y(^2) " 0y(^2)) (30J'(X3) + 0j(x3))] 

■Eg{tg)hg{A, B, C, b^, ^2 , ^3 , X 2) " Eg ( t'^ ) /i, ( A' , B\ C\ ^2 , ^1, ^3 , X {) 

■ [-4rv0^(a;2X(x3) + 2rvr,X30t'(x2) (3<(a;3) + <(a;3))] } , (46) 
for the case of -B ^ V" transition, and 

Ml^y = —=GFCpm% / dxidx2dx3 / bidbib2db2(pB{xi,bi) {[- {1 - X2) {2(t)^{x2) + Tn 
v3 Jo Jo 

■(30J'(X2) - 0J'(X2)) + r^X2 (0^(X2) + (I)J,{X2)) } ^^(Xg) + Ty (1 + X2) X30;^(X2) 

■ (30^(X3) + 4>v{^3)) + r^rv (1 - X2) (0j'(x2) + <l)^{x2)) {^M^s) - (pvixs)) 
+rr,rvX3 (1 - 2x2) «(X2) - 0^(0:2)) (30t-(S3) + M^s))] 

■Eg{tg)hg{A, B, C, b^, ^2 , ^3 , X 2) " E g {t'^) h g {A' , B', C', ^2 , ^1, ^3 , ^^ 1 ) 

■ [4r^^^{x2)(pv{x3) - 2r^rvX3<(x2) (30t'(x3) + ^^3))] } ■ (47) 

for the case oi B ^ rj transition. Here the hard scale tg and t' are the same as in Eq. (HTl) . 
The evolution factor Eg{t) in Eqs. fH^ and fH7|) is of the form 

^9(^) = C|f (i) «'(^) ■ exp [Smg] , (48) 
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with the Sudakov factor Smg and the hard function hg, 



Smgif) = s \ximB/V2, bA + s \X2mB/V2, 62J + s n 1 - X2)mB/V2 
+s (x3mB/V2, hj + s Ul- xs)mB/V2, h. 



^3 
1 " 



A 



j^ ln(t/A) ^ j^ ln(t/A) ^ ^^ ln(t/A) 



-ln(6iA) -ln(62A) -ln(63A)_ 



(49) 



hg{A, B, C, h, h, h, X,) = -St{x,) Ko{Bh) KoiCh) 

'.7r/2 

d6'tan6' ■ Jo(^&i tan6') Jo(^&2 tan6') Jo(^&3 tan6'), (50) 



where the functions -^'0(2;) and ^0(2;) are the Bessel functions, the form factor St{xi) with 
i = 1,2 has been given in Eq. flB7p . and the invariant masses A^'\ B^'^ and C*^'^ of the 
virtual quarks and gluons are of the form 



A = y/x^rriB, B = B' = y/x^niB, C = 2a/(1 - X2) 



XsTUb, 



A' = y/xlniB, C" = v|xi -xslms. (51) 

For B — > Vrj' decays, we find the similar results by making appropriate replacements. 

The total "chromo-magnetic penguin" contribution to the considered B — * Vr]^'^ decays 
can therefore be written as 



'y,7(') ~ ' 'I \"-eff\ 



<:? = < Vr^^'^\n-^\B >= A, M^^) + M^Jl , (52) 



(9) , l\A9) 



Vri(') ^ riOV 



where At = VtbV^*^. Again, the chromo-magnetic penguins do not contribute to i? ^ (pi]^'' 
decays. 

V. NUMERICAL RESULTS AND DISCUSSIONS 

Using the wave functions and the central values of relevant input parameters as given 
in Appendix El we firstly find the numerical values of the corresponding form factors at 
zero momentum transfer: 

A^-^iq' = 0) = 0.32tlfM), 
<-(g^ = 0) = 0.29+°:°^(a;,), 
^s^^(')^^2 = 0) = 0.22 ± Om{uJb), (53) 

for ujh = 0.40 ± 0.04GeV, which agree well with those obtained in QCD sum rule calcula- 
tions. 
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A. Branching ratios 

For a general charmless two-body decays B — >■ Vr]^'\ the branching ratio can be written 
in general as 



Br(B -^ Vr] 



W^ 



Tb 



m 



IQuniB 
where tb is the lifetime of the B meson, and the decay amplitude is the form of 



M =< Vv^'^lHeff + Kff + 7^™ 15 > . 



Acmp) I 



(54) 



(55) 



Using the wave functions and the input parameters as specified in previous sections, it is 
straightforward to calculate the CP-averaged branching ratios for the considered decays, 
which are listed in Table [B For comparison, we also list the corresponding updated 
experimental results [22, l23| an d numerical results evaluated in the framework of the 
QCD factorization (QCDF) 



231 an c 



TABLE I: The pQCD predictions for the branching ratios (in unit of 10^^). The label LOnlowc 
means the LO results with the NLO Wilson coefficients, and +VC, +QL, +MP, NLO means the 
inclusion of the vertex corrections, the quark loops, the magnetic penguin, and all the considered 
NLO corrections, respectively. 



Mode 


LO 


LOnlowc +VC +QL +MP 


NLO 


Data 


QCDF 


i?±- 


* P'^V 


6.9 


7.4 


6.8 


7.5 


7.2 


6.7 


5.4 ±1.2 


9.41^:^ 


B^- 


>pV 


5.2 


4.8 


4.6 


4.9 


4.7 


4.6 


q 1+3-7 
^•^-2.8 


6.3it^ 


B^^ 


p^r, 


0.08 


0.08 


0.19 


0.16 


0.12 


0.13 


< 1.5 


0.03l°i^ 


B^^ 


pS 


0.05 


0.04 


0.13 


0.06 


0.04 


0.10 


<1.3 


r,r,i +0.12 


B^^ 


Lvr] 


0.22 


0.34 


0.67 


0.33 


0.25 


0.71 


<1.9 


n qi+o-46 


B^^ 


ujr]' 


0.12 


0.18 


0.52 


0.19 


0.15 


0.55 


<2.2 


20+°-^*^ 


B^^ 


<Pv 


0.001 


0.002 


0.011 


- 


- 


0.011 


<0.6 


0.001 


B^^ 


H 


0.096 


0.053 


0.017 


- 


- 


0.017 


<0.5 


0.001 



It is worth stressing that the theoretical predictions in the pQCD approach still have 
relatively large theoretical errors induced by the large uncertainties of many input pa- 
rameters, such as a;b, Gegenbauer coefficient 02, the CKM angle a and m^. The pQCD 
predictions with the major theoretical errors for the branching ratios of the decays under 
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FIG. 5: The a dependence of the branching ratios of B'^ -^ p'^rj^'' decays for ujb = 0.36 GeV 
(dotted curve), 0.40 GeV (solid curve) and 0.44 GeV (dashed curve). 



consideration are the following 



Bri B^ - 

Br{ B^ - 

Br{ B° - 

Br{ fi° - 

Br{ fi° 

Br{ 5° - 

Br{ 5° 

Bri B^ - 



P^V') 
-A) 
■ P'V') 



[6.7+?:^(o.,);i;°(/.o)+Si(«)l^i(a2)+^:^,(a2p)] x 10 



-6 



1.4 

1.1 

4 

2 

3 

2 

1.7 
1.3 

1.3 
1.1 



' , ^+o.5 

.^feJ-0.7 



{fio)t?,iia] 



-0.3V 



\+0.1 



.^b)-OAi^^o)-0.o{'^)-0.o{(^'2)Lo.l{0'2p 



Ub 



-oaVP'O, 



±0.2(«^+oo 



[4.6 
[1.3+°: 

[1.0+°: 

[7.1 
[5.5 

[i.i±o.i(^,)+^:^(/.o) 

[l.7±0.2{uJbnifMt'oi{^s 



±0.4(a2) 



l-0.0/_ 
-0.ll^2p 



\+0.2/ 



X lo^^ 

X 10"^ 
^5:j^a2j±0.1(a2p)] X 10' 



\+l-2/ 



+'-'(«2)+^:^(a2.)] X 10 



-7 



-1.41 



+0.1 



(^.)^L6(/^o)^l:^(«)^l:^(a2)l[]:?(a2.) x lo 



1-7 



+0.4, 
-0.2> 



rrio 



hO.l 
-0.2 



1-8 



± 0.1(02) ±O.l(a20)] X 10~^(56) 



where the major errors are induced by the uncertainties oiuf, = 0.4±0.04 GeV, /xq = 1.0± 
0.5 GeV, a = 100° ±20°, m, = 130 ±30 MeV, Gegenbauer coefficients 02 = 0.115 ±0.115, 
ct2p = Ci2ui = 0.15 ± 0.15 and a2(j> = 0.2 ± 0.2, respectively. 

In Figs. [5l M and [7] we show the a and tUfe-dependence of the pQCD predictions for 
the branching ratios of -B ^ pil^'\ B -^ lut]^'^ decays for cub = 0.4 ± 0.04 GeV, and 
a = [0°, 180°], a2 = 0.115 and a2p = a2^ = 0.15. 

From the numerical results and the figures, we observe that 

• For B^ — ^ p^rj decay, the inclusion of the considered NLO corrections can improve 
the agreement between the pQCD prediction and the data. But for B^ — > p^r/' 
decay, we are not so lucky. Although the pQCD predictions for Br(B± -^ p^r]^'^) 
agree with the data within one standard deviation, but the predicted pattern of 
Br{B± — »• p^rj) > Br{B^ri') in both the pQCD and QCDF is contrary to the 
observed one. 



For i?° — »• p^ [oj ^ (f))ri^'' decays, the pQCD predictions for their Br's are consistent 
with currently available upper limits. Except for Br{B -^ (prj'), the inclusion of the 
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FIG. 6: The a dependence of the branching ratios (in units of 10 ^) of B^ — > p^Tj^'' decays for 
LOh = 0.36 GeV (dotted curve), 0.40 GeV (sohd curve) and 0.44 GeV (dashed curve). 
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FIG. 7: The a dependence of the branching ratios (in units of 10 ^) of B^ -^ mr]^'' decays for 
ujh = 0.36 GeV (dotted curve), 0.40 GeV (sohd curve) and 0.44 GeV (dashed curve). 



partial NLO contributions to other decays can enhance their Br's by a factor of two 
to ten, and generally larger than the QCDF predictions, which will be tested by the 
forthcoming LHCb experiment. 

B. CP-violating asymmetries 

Now we turn to the evaluations of the CP-violating asymmetries of i? — > p{uj, (f))r]^'^ 
decays in pQCD approach. For B~^ -^ p^rj and B'^ — > p+ry' decays, the direct CP-violating 
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asymmetries Acp can be defined as: 

The pQCD predictions for tlie direct CP-violating asymmetries of B"^ -^ p^^i]^'^ decays 
are listed in Table HTl For comparison, we also list currently available experimental results 



22 . |23| and the numerical results evaluated in the framework of the QCD factorization 



(QCDF) [19|. 

TABLE II: The pQCD predictions for the direct CP-violating asymmetries of B -^ p t]^'' 
decays (in units of 10"^). 



Mode 


LO +VC +QL +MP 


NLO Data QCDF 




- P^V') 


0.0 1.3 1.4 -0.1 
-6.8 -25.3 -5.7 -7.1 


1.9 1.0 ±16.0 2.4 
-25.0 -4.0 ±28 -4.1 



The NLO pQCD predictions for the central values of the direct CP-violating asymme- 
tries and the major theoretical errors for 5^ -^ p'^rf^'^ decays are 

A'^},{B^ - P^V) = [l.9lg:^K)lgJ(«)ig:^(a2)+g|(a2,)] x 10-^ (58) 

A'^liB^ - P^V) = [-25.0+°:^(.;,)+li(a)+°;?(a2)+?i(a2,)] x 10-^ (59) 

where the major theoretical errors come from the variations of Uf, = 0.4 ± 0.04 GeV, 
a = 100°±20°, Gegenbauer coefficients oa = 0.115±0.115, agp = a2uj = 0.15±0.15. Both 
the pQCD and QCDF predictions are consistent with the data because of the still large 
theoretical and experimental errors. In Fig. [8], one shows the a and ti^fe- dependence of the 
LO and NLO pQCD predictions for the CP-violating asymmetries of S^ -^ p^ri^'\ 

As to the CP-violating asymmetries for the neutral decays i?° -^ p'^(uj)'r]^'\ the effects 
of B^ — B^ mixing should be considered. The CP-violating asymmetries for such decays 
are time dependent and can be defined as 

P (50(At) ^ fcp) - P (5S(At) -. fcp) 

P (50(At) -. fcp) + P (i?°(At) -. fcp) 
= Afp cos(AmAt) ± A™)? sin(AmAt), (60) 

where the direct and mixing induced CP-violating asymmetries A'^'^p and A^'^p can be 
written as 

Adir _ \^CP\ - 1 ^mix _ "^ImjXcp) , . 

-L ± \^CP\ -L ± \^CP\ 

with the CP-violating parameter Xcp is 

_ V,lVu{Vv^'HH,ff\^) 

'^' - -v,,v:,{vvmeff\Br ^''^ 
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FIG. 8: The a an a;f,-dependence of the CP-violating asymmetries of B — > p r]^'' decays for 
a = [0°, 180°] and Ub = 0.36 GeV (dotted curve), 0.40 GeV (sohd curve) and 0.44 GeV (dashed 
curve) . 



TABLE III: The pQCD predictions for the direct, mixing-induced and total CP asymmertries 
of B° -^ pO(^)^(/) decays (in unit of lO^^). 



Mode 


LO -LOnlowc 


+VC 


+QL +MP 


NLO 


QCDF 


A^SUB' - 


-.A) 


-78.3 


-94.4 


-91.4 


-79.0 


-97.4 


-89.6 


— 


Ai^UB' - 


^AO 


77.3 


-42.6 


-79.8 


-96.8 


-86.2 


-75.7 


- 


A'SUB' - 


■^ijjr\) 


94.6 


53.4 


45.7 


48.7 


32.1 


33.5 


33.4 


A'^'piB^ - 


^ur]') 


30.0 


26.9 


26.5 


15.4 


-12.1 


16.0 


-0.2 


A'^fiB^ - 


-^A) 


44.7 


26.1 


16.9 


14.1 


17.5 


22.7 


- 


A^}^{B^ - 


-. p^) 


-24.0 


-15.9 


-46.0 


-3.5 


-14.4 


-49.0 


- 


A^']?{B^ - 


-^ UT]) 


-7.4 


83.8 


40.3 


81.5 


80.2 


39.0 


- 


A'^^iB^ - 


-^ui) 


58.9 


88.7 


78.6 


82.2 


71.6 


77.0 


- 


^*?M^° - 


^A) 


-24.0 


-46.4 


-48.9 


-42.5 


-52.4 


-45.0 


- 


A'S'piB'' - 


. pS) 


35.5 


-34.3 


-72.1 


-62.2 


-60.8 


-71.0 


- 


^*?M^° - 


^LUT]) 


54.2 


73.9 


48.1 


69.9 


58.9 


39.8 


- 


A'S'p^B^ - 


■^ui) 


48.4 


59.7 


54.6 


49.4 


27.1 


47.3 


- 



If we integrate the time variable t, we will get the total CP asymmetries for B^ — > Vrj^'^ 
decays, 



A 



CP 



AfT^ + 



X 






l+3;2-CP-^^^2-CP, 



(63) 



where x = Am/F = 0.775 for the B^ — B mixing 22 . 

The pQCD predictions for the CP-violating asymmetries and the total CP violation 
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when different NLO contributions are included step by step are listed in Table III. The 
pQCD predictions with major theoretical errors are given in Eqs. 






-2 

-2 



^^]?(5°^A) = [22.7±6.lK)™(«)+?2'5(a2)li:^(a2p)]xl0- 

^—(5°^ A') = [-49.0lJ|K)ir(a)+ii(a2)+l?:^(a2,)] X 10- 

^^-(5° - cov) = [39.0+°i(.;,)+^°:^(a)+i^(a2)+?:^(a2j] x 10-^ 

A^'p^iB' ^ 0.V) = [77.0+°:f(a;,)+i:°(a)+°:?(a2)l°:^(a2.)] x 10-^ (65) 



3-^ 

^-2 

5 



■A*c°M5°-A') = [-7i-o+l?K)l^:^(«)^l:?(a2)l^:^(a2,)] X ir^ 

^- (i?° - c.r^) = [39.8;°;e^(a;,);^?J(a)l|^(a,)+|^(a,J] x 10-^ 

^*e«M5° - ^V') = [^7-^tUi^b)lliAc^)tliia2) ± 1.2(a2.)] x 10-^ (66) 

where the dominant errors come from the variations oiuf, = 0.4±0.04 GeV, a = 100°±20°, 
and Gegenbauer coefficient a2 = 0.115 ± 0.115, a2p = a2w = 0.15 ± 0.15. 

For the CP- violating asymmetries of B^ — >■ p^(uj)'r]^'^ decays, unfortunately, there is no 
data available currently. For B^ -^ (pr]^'^ decays, there is no CP violation. The reasons are 
simple: (a) the total decay amplitude at the LO level as given in Eq. flTTl) is proportional 
to only one CKM factor C,t', and (b) among the NLO contributions considered here, only 
the vertex correction ( real correction ) is relevant for this decay mode. 

VI. SUMMARY 

In this paper, we calculate some NLO contributions to the branching ratios and CP- 
violating asymmetries of 5^ -^ p^^i]^'^ and B^ -^ p^{uj,(f))r]^'^ decays by employing the 
pQCD factorization approach. 

From our calculations and phenomenological analysis, we found the following results: 

• The pQCD predictions for the form factors oi B ^ p,uj and i]^'^ transitions are 

A^-^^(O) = 0.32+°:°^(a;b), <^"(0) = 0.29+°;°^ (a;,,) and F^-'''"\o) = 0.22±0.03(a;b) 
for LJb = 0.40 ± 0.04GeV, which agree very well with those obtained in QCD sum 
rule calculations. 

• For B^ —>■ p^T] decay, the inclusion of partial NLO contributions can improve the 
agreement between the pQCD predictions and the measured values. For the neutral 
decays, the NLO contributions can provide significant enhancements to the LO 
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predictions: 

Br{ B^ -^ p^r]) = [6.7+?;^] x 10"^ 
Br{ B^ -* p^r]') = [4.6+li] x 10~^ 

5r(i?°^A) = [1.3+^x10-^ 

Br{ B^ -^ p^T]') = [1.0 ± 0.5] X 10"^ 

Br{B'^ur]) = [7.ltli] x 10"^ 

Br{ 5° -^ LUT]') = [5.5tli] x 10"^ 

Br{ 5° -^ (j)r]) = [l.l+ag] x 10~^ 

Br{ 5° ^ 0r/') = [l.7t}V] x 10-^ (67) 

where the various errors as given in Eq. (l56|l have been added in quadrature. 

• The pQCD predictions for AfipiB"^ — > p'^i]^'^) are consistent with the data, but both 
the theoretical and experimental errors are still large. For other neutral decays, the 
pQCD predictions for CP violating asymmetries are generally large in magnitude 
and could be tested by the forthcoming LHCb experiments. 

• Only the NLO contributions from vertex correction, quark-loops and chromo- 
magnetic penguins are calculated here. The NLO corrections from the hard- 
spectator and annihilations diagrams are still absent now. It is an urgent task 
to do the relevant calculations, in order to provide a complete NLO calculation in 
the pQCD approach. 
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APPENDIX A: WAVE FUNCTIONS AND INPUT PARAMETERS 

The B meson is treated as a heavy-light system. For the B meson wave function, since 



the contribution of (p^ is numerically small 24] , we here only consider the contribution 
of Lorentz structure 

\/ZI\c 



with 



(j)B{x,b) = Nbx'^{1 — x^exp 



Ml x^ I, ,,; 

— V - oM) 



2uji 2 



(A2) 



where ujb is a free parameter and we take Ub = 0.4 ± 0.04 GeV in numerical calculations, 
and Nb = 101.445 is the normalization factor for Uf, = 0.4. 

21 



For the considered decays, the vector meson V is longitudinally polarized. The longi- 
tudinal polarized component of the wave function is defined as: 



V^K 



{/ [mv(l)vix) +^v(f)vix)] +mv0y(2;)} 



(A3) 



where the first term is the leading twist (twist-2) wave function, while the second and 
third terms are twist-3 wave functions. 

The twist-2 DA's for longitudinally polarized vector meson can be parameterized as: 



(pv{x) 



h 



2^/2N, 



:6x(l-x) l + a2vC^2 {^x-l) 



3/2, 



(A4) 



for V = p,uj,(f)] and fy is the decay constant of the vector meson with longitudinal 
polarization, and numerically [25|: 



L = 216MeV, f^ = 187MeV, 



215MeV. 



(A5) 



The Gegenbauer coefficients have been studied extensively in the literature. Here we 
adopt the following values from the recent updates |25|]: 



a2p = a2<^ = 0.15 ±0.15, a2<^ = 0.2 ± 0.2. 



(A6) 



We shall vary the Gegenbauer coefficients of the twist-2 distribution amplituds by 100%, 
which is larger than the error specified in 25||. Therefore, the theoretical uncertainty of 
our predictions from this source is conservative. 



As for the twist-3 DAs (py 
0y(x) 



and (/)*/, we adopt their asymptotic form 26|: 



^^^ (l-2a;), <p'y{x) 



^^^ -.{2x 



For 77*^') meson, the wave function for qq (g = -u, d) components of rf^'^ meson are given 



as 



%^{P,xX) = 



V2N, 



=75 



c 



Vq' 



x) +ml''^^(x) +C"^o' 



v-n)(l)'^(x) 



(Ai 



where P and x are the momentum and the momentum fraction of rjq, respectively. We 
assumed here that the wave function of rjq is same as the vr wave function. The parameter 
( is either +1 or —1 depending on the assignment of the momentum fraction x. The ss 
component of the wave function can be similarly defined. 

For the mixing oi rj — rj' system, we here use the quark-fiavor basis, that is the T]q = 
{uu + dd)/ \/2 and rfg = ss. Then the physical states r] and i]' are related to the fiavor 
states through a single mixing angle (j), 



cos (f) — sin I 
sin 6 cos (h 



Vg 

Vs 



(A9) 



The relation between the decay constants (/^, /^, /^,, /^,) and {fg, fg, ) can be written as 

f^ = fqCos(j), /^ = -/,sin0, 

f^, = /,sin0, /^^, = /,cos0. (AlO) 
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The chiral enhancement ttIq and m'^ associated with the two-parton twist-3 1]^ and rjs 
meson distribution amphtudes have been defined as |15i] 



m 



rrin 



rrin 



qq 



2mq 

mis 
2r7io 



\m cos (p + m„/ sm — 

zrrin 



- — [ml, cos^ (p + ml sin^ < 
2m^ ' ' 



fq 



^ 'ml^i — m^) cos (j) sin i 



fs 



'' 'm^, — m^) cos sin < 



(All) 
(A12) 



by assuming the exact isospin symmetry rriq = rriu = rrid- The three input parameters 
fq, fs and (f) have been extracted from the data of the relevant exclusive processes 271]: 

/, = (1.07 ±0.02)/,, /, = (1.34 ±0.06)/,, = 39.3° ± 1.0°, (A13) 

It is still unclear for the possible gluonic component of 7]' meson. From currently known 
studies I id. I 111 . [l2] we believe that there is no large room left for the contribution due to 
the gluonic component of?]', and therefore will neglect the possible gluonic component in 
r)' mson. 

The distribution amplitude 4>^'^'^ represents the axial vector, pseudoscalar and tensor 
component of the wave function respectively |23] • They are given as: 



<(^) 



(PvM) 



(P'^S^) 



^"^ 6x(l - x) [l ± afCl^\2x - 1) ± a^'Cf '(2x - 1) 



2y/2N, 

+al''Cl^\2x - 1 

J-Hq 



2^/2N, 
-3 

JVq 



I + {SOvs - -pi)C',^\2x - 1) 
3 |^3C^3 ± 7^pL(1 ± ea^)) Cy\2x - 1) 



2Q^vq 

, ;i-2x) 

■ (l-10x±10x2)] 



1 7 3 

1 ± 6 ( 5r73 - -r/3^3 - —p\ - -p\/2 



with 



"qqi ^\ 



Vqq 



Pn^ = 2mq/m 
and the Gegenbauer polynomials C^(t) 

'2 
-y3/2 

■y3/2 



0, a"","" 



0.115 ±0.115, al" 



-0.015. 



crit) = -{3t'-i) 



■3/2/ 



cnt) = 3t, crw = ^(5t'-i) 



C]/'(t) = l(3-30t2±35t^), 

o 

3 



Cl'^t) 
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i-i4r ±2ir). 



(AM) 

(A15) 

(A16) 
(A17) 

(A18) 
(A19) 
(A20) 



The Gegenbauer coefficients can vary by 100%, but we do not consider the uncertainty 
from the coefficients al'''' , to which our predictions are insensitive. The values of other 
parameters are 773 = 0.015 and uj = —3.0. 
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As to the wave function of the ss components, we also use the same form as qq but 
with some parameters changed : 

Prj, = 2ms/mss, a?^ = 0, a^= = 0.115 ± 0.115, aj' = -0.015. (A21) 

Besides those specified in the text, the following input parameters will also be used in 
the numerical calculations: 

/^ = 130MeV, fn = 210MeV, m^ = 547.5MeV, m^, = 957.8MeV, 

rUg = 5.6MeV, m, = 130 ± SOMeV, uib = 5.28GeV, 

rup = 774MeV, m^ = 780MeV m^ = 1.02GeV, mw = 80.41GeV, 

tbo = 1.528ps, tb+ = 1.643ps, (A22) 

For the CKM quark-mixing matrix, we use the Wolfenstein parametrization as given 



in Ref. [22,123 



Kd = 0.9745, Vus = X = 0.2200, |Kb| = 4.31 x 10~^ 
Vcd = -0.224, Vcd = 0.996, V^b = 0.0413, 
\Vtd\ = 7.4 X 10-3, Vts = -0.042, \Vtt\ = 0.9991, (A23) 

with the CKM angles p = 21.6°, 7 = 60° ± 20° and a = 100° ± 20°. 

APPENDIX B: RELATED FUNCTIONS 

We show here the function /ij's, coming from the Fourier transformations of H^^\ 

he{Xi,X2, bi, 62) = Kq {y/x^niBbl) [9{bi - b2)Ko (y/x^niBbl) h (v^"^iJ&2) 

+9{b2 - bi)Ko {^/x^mBb2) h (v^^^s^i)] St{x2), (Bl) 

ha{x2, X3, 62, ^3) = Kq [i\/ll^^X2Jx^mBb2J [6{b3 - 62)^0 {i^/x^mBb^) Iq (« Va^"^i?^2) 

+9{b2 - 63)^0 {iy/x^mBb2) Iq (iy/x^niBbs)] St^xs), (B2) 



/l„(Xi,X2, 0:3,61,63) = <^ e{bi - b3)Ko{MBy/XiX2bi)lo{MBy/XiX2b3) 
+ 9{b3- bi)Ko{MBy/XiX2b3)loiMBy/XiX2b^ 



yHo(a/(x2(x3 - xi))MBb3), for xi - X3 < ■ , . 
K}^\^/ix^(x^'^^xs)MBbs), for xi - X3 > ' 



hnaixi,X2,X3,bi,b3) = \ 9{bi - b3)Ko{i^/{l - X2)x3biMB)loii^/il - a;2)a;363Mj 



B 



+ (^(63 - 6i)Ko(V(l - X2)x363Mb)Io(V(1 - X2)x3biM, 



B, 



Kq{Mb\/{xi - X3)(l - X2)6i), for xi - X3 > , ,^^. 
^'-}l^^\MB^{x3-x{){l-X2)b^). for Xi - X3 < 0/ '^ ^ 
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+ e{b3 - 6i)Ko(2 V(l - X2)x3hMBMty/{l - X2)x3biMB) 

Ko(MsFA), for F^ > o' 
f}l'i\MB^AW\bi), for F2<0^ 



(B5) 



where Jq is the Bessel function and Kq, Iq are modified Bessel functions Kq^—ix) 
— {n/2)YQ{x) + z (77/2) Jo (a;), and i^(i)'s are defined by 



"h 



1 -X2(l -X3 -Xi). 



(B6) 



The threshold resummation form factor St{xi) is adopted from Ref .[28!]. It has been 
parametrized as 



2i+2T(3/2 + c) 



(B7) 



where the parameter c = 0.3. This function is normahzed to unity. The evolution factors 
Ee and Ea are given by 

Ee{t) = as(t)exp[-Sab{t)], 

Kit) = a,(t)exp[-5erf(t)]|62=bi, 

Ea(t) = as(t)exp[-Sgh{t)], 

E'^it) = a,(t)exp[-5e/(t)]|b,=63, (B8) 

The Sudakov factors used in the text are defined as 

Sab{t) = s (ximB/V2, bij + s (x2mB/V2, 62J + s n 1 - X2)mB/V2, h 



1 

"a 



In 



ln(t/A) 



In 



ln(t/A) 



-ln(6iA) -ln(62A) 

Scd{t) = s (ximB/V2, bi) + s (x2mB/V2, bA + s Ul - X2)mB/V2, &i 

+s (x3mB/V2, 63) + s Ml - X3)mB/V2, k 
1 



(B9) 



Pi 



21nJ^^(^ + ln.l^(^/^) 



ln(6iA) -ln(63A) 

Sefit) = s \^XimB/V2,bi] + s (x2TnB/V2,b2] + s Ul- X2)mB/V2,br 

+s (x3mB/V2, 62J + s n 1 - X3)mB/V2, 62] 

_J_r ln(t/A) ln(t/A) 

A [ -ln(6iA)^ -ln(62A) 

Sgh{t) = s (x2mB/V2, 62] + s (x3mB/V2, b 
+ s (^(1 - X3)mB/y/2, 63 j - — 



(BIO) 



s nl -X2)mB/V2,k 



l^J^(VAl_^l^,ln(t/A) 



ln(6iA) 



ln(62A) 



(Bll) 



(B12) 
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where the function s{q, b) are defined in the Appendix A of Ref.p]. The scale tj's in the 
above equations are chosen as 



ta 


= 


niax( 


t'a 


= 


niax( 


tb 


= 


niax( 


ti 


= 


niax( 


tc 


= 


max( 


t'c 


= 


max( 


td 


= 


niax( 


t'd 


= 


niax( 



X2mB, y/XiX2mB, I/61, I/62) , 

xlniB, ^/x^ruB, 1/&1, 1/&2) , 
a/x2|1 -X3 -xilms, ^Jx^ruB, l/h, I/63) , 
\/x2\x's, -Xi\mB, y/x^niB, 1/61, 1/63) , 
a/(1 -X2)x3mB, a/|xi -x-i\{l -X2)mB, 1/61, 1/&3) , 
a/|1 - 2:2(1 - X3 - xi)\mB, a/(1 - X2)x3mB, I/61, 1/&3) , 

a/(1 -X2)x3mB, a/(1 -X2)mB, 1/&2, I/&3) , 

v/(l -X2)x3mB, y/x^ms, 1/&2, I/&3) • (B13) 
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